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Abstract 

Electric dipole moments (edms) of bound states that arise from the constituents having edms 
are studied with field-theoretic techniques. The systems treated are the neutron and a set of 
paramagnetic atoms. In the latter case it is well known that the atomic edm differs greatly 
from the electron edm when the internal electric fields of the atom are taken into account. In 
the nonrelativistic limit these fields lead to a complete suppression, but for heavy atoms large 
enhancement factors are present. A general bound-state field theory approach applicable to both 
the neutron and paramagnetic atoms is set up. It is applied first to the neutron, treating the quarks 
as moving freely in a confining spherical well. It is shown that the effect of internal electric fields 
is small in this case. The atomic problem is then revisited using field-theory techniques in place 
of the usual Hamiltonian methods, and the atomic enhancement factor is shown to be consistent 
with previous calculations. Possible application of bound-state techniques to other sources of the 
neutron edm is discussed. 
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I. INTRODUCTION 



One of the strongest constraints on new physics beyond the standard model is provided 
by measurements attempting to detect a nonvanishing electric dipole moment (edm) of the 
neutron, heavy atoms, or molecules This is because the extremely small value of the edm 
of quarks and the electron in the standard model means that any detection with presently 
available sensitivities arises from such new physics. According to Schiff's theorem |2J, the 
edm of atoms or molecules arising from an electron edm vanishes in the nonrelativistic limit. 
However, the power of experiments on heavy atoms or molecules to put limits on a possible 
electron edm d e is in fact increased by orders of magnitude because of a large violation of 
Schiff's theorem discovered by Sandars which stems from relativistic effects. Typical 
values of the enhancement factor R, defined by 

^atom = Rd e: (1) 

are R = 120 for cesium 4] and R = —685 for thallium These calculations of R are 
carried out in a Hamiltonian formalism, with the relativistic many-body problem treated 
using many-body perturbation theory (MBPT), generally summing infinite classes of MBPT 
diagrams for higher accuracy. 

One of the purposes of the present paper is to recalculate the enhancement factor for 
paramagnetic atoms, specifically thallium and the alkali metals lithium through francium, 
using a different formalism, one which treats corrections in terms of one photon exchange. 
Our aim here is not to achieve high accuracy in the many-body part of the problem, but 
rather to formulate the calculation in a field-theoretic framework |6| . This framework allows 
us to address the main purpose of the paper, calculating the neutron edm induced when 
the up and down quarks have edms d u and dd respectively. We are specifically concerned 
with the question of how the internal electric fields arising from the charged quarks affect 
the edm of the bound state, since their role is so important in the atomic case. In order 
to apply the techniques of atomic physics to the neutron, we treat the up and down quarks 
as moving freely in a confining well. This is a simplified version of the MIT bag model 
[t| , without complicating issues such as bag pressure and the nonlinear boundary condition; 
we also neglect gluon exchange. We refer to this approach as the static well model in the 
following. The model is known to give the neutron (and proton) magnetic moment to within 
20 percent of the experimental value when the quarks are treated as nearly massless particles 



with magnetic moments coming from their confinement in the well [8|, so using it to evaluate 
the neutron edm that arises from nonvanishing values of d u and dd should give results of 
similar accuracy. 

The paper is organized as follows. In the next section we introduce the bound-state field 
theory tools that will be used for both the neutron and paramagnetic atom calculations, both 
of which are variations of the Furry representation [9]. In the following section we calculate 
first the lowest-order edm of the neutron, finding a 17 percent reduction from the standard 
nonrelativistic result, and then evaluate the effect of one photon exchange between quarks, 
which is equivalent to calculating the effect of the internal electric fields in the neutron. This 
effect will be shown to be quite small, so that despite the highly relativistic nature of the 
quarks within the model, we do not find large edm enhancement factors for the neutron. 
We then turn to the atomic problem, beginning with a brief description of the violation 
of Schiff's theorem [sj in atomic physics for paramagnetic atoms, after which we present a 
field-theoretic treatment of the problem analogous to that used for the neutron. We show 
that the approach leads to a suppression for light alkalis along with the well-known large 
enhancement factors for heavy atoms. In the concluding section we discuss how different 
sources of CP nonconservation could be treated with the methods described here, and also 
ways in which the atomic calculations could be improved. 

Because we will be dealing with elementary particles with different electric charges, in 
the following we will present all formulas with factors of q a , with q e = — e, q u = 2/3e, and 
qd = — l/3e, e taken to be positive. 

II. BOUND-STATE FIELD THEORY 

For both problems treated here we assume spherical symmetry, and we expand field 
operators in terms of solutions to the Dirac equation 



with energy eigenvalues e n , and where Xk/i is a spherical spinor. To treat the bound state 
problem we use variants of the Furry representation, originally introduced to evaluate QED 
effects in hydrogen. This representation can be thought of as intermediate between the 
standard interaction and the Schrodinger representations. 
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To transform from the Schrodinger to the interaction representation the unitary trans- 
formation 

\Mt)} = e^lMt)} (3) 

is carried out, with 

K-JWVWt+finW.*)- (4) 
In the Furry representation, the unitary transformation is 

\Mt)) = e^lMt)}- (5) 

where, defining r = \x\, 

/Z oi 
d 3 ''xij)*(x)[d ■ p '+ 0m(r) (6) 

We have replaced the standard electron mass term j3m with a position-dependent mass 
for later generalization to the static well model, though it is of course constant in atomic 
physics. This representation is used to this day to treat radiative corrections in hydrogen 
and hydrogenic ions, where the primary role of the nucleus is to provide a classical Coulomb 
field centered at the origin, which introduces the extra term into H . (Finite nuclear size 
effects, which are particularly important for highly charged ions, are accounted for by taking 
Z — > Z(r), with Z(r) modeled by a Fermi distribution). 

After the transformation to the Furry representation the interaction Hamiltonian (for the 
atomic problem) has the usual form 



Hj = q e j d 6 xiP(x)^(x)A»(x) (7) 

with N 1 the quantized radiation field, but with the electron field operators now expanded 
in terms of solutions to the Dirac equation in an external Coulomb field of the form given 
inEq. ^. 

To evaluate e nerg y corrections, a generalization of the Gell- Mann-Low approach intro- 



duced by Sucher [10j] can be used 



AE= lim -A ln (0l r (e" iA/d:Eoe " E|a:olH/(:Eo) )l0), (8) 

<e->0,A->1 2 OA 

which has an S'-matrix modified by the factor e~ e ' x °', included in order to adiabatically 
turn off the interaction at large positive and negative times. This allows standard Feynman 



diagram techniques to be applied, with the adiabatic factor usually leading trivially to a 
factor 1/e that cancels the e in the numerator of the above formula. 

In this work we are interested in treating many-electron atoms and the neutron, but only 
slight modifications of the approach used for hydrogen described above are needed. We start 
with the neutron. We model it as three light quarks confined in a spherical well of radius 
R = 1.2 fm, choosing the scalar mass of the quarks to be finite inside the well and infinite 



outside, which leads to the boundary conditions of the MIT bag model [7]. Just as with 
the description of hydrogen, we note that there is a fixed origin, in this case the center of 
the well instead of the position of the nucleus. In both cases this is equivalent to neglecting 
recoil effects. Because the masses of the up and down quarks we will be considering are 
very light, we set them equal to zero. The calculations presented below were tested with 
nonzero quark masses and found to be insensitive to this approximation. A simplification of 
choosing the up and down quark masses to be equal is that the spatial wave functions are 
then identical. 

For the neutron, then, we take 

H = Xw d 3 x[ip ui \x)[a ■p + r rom(r)]i/j U i(x) + ipd^ (x)[a ■ p + "/ m(r)]ip di (x)] (9) 

i 

where i is a color index. In this case the quark fields are expanded in terms of solutions 
of the free Dirac equation for massless particles with MIT bag model boundary conditions. 
The neutron wave function consists of a product of three Is wave functions [8], with 

= ^k[-^\ b b) c bl c + ij>)Ac + blb] c blM, (10) 



2VU 

where we have introduced a notation in which a and b denote spin up and down states of an 
up quark, and c and d spin up and down states of a down quark. The color indices and 
k are again understood to be summed over. The ground state energy of a massless quark 
when R = 1.2 fm is 335.9 MeV, which will be denoted e g in the following: we also denote 
the associated wave function as ip g {x). This defines the starting point for our treatment of 
the neutron edm, and calculations involving the appropriate interaction Hamiltonians can 
now be carried out in a systematic manner. 

We now turn to the atomic case, where we will be interested in treating atoms with many 
electrons, in which case the original Furry representation is inadequate. To approximate the 
effects of screening, we introduce a local, central potential chosen to be close to the Hartree- 
Fock (HF) potential. Most atomic calculations based on Hamiltonian methods in fact use 



the latter potential, but its nonlocality makes it unsuitable for field theory calculations. The 
new free Hamiltonian is 

/Z ot 
d 3 xip Jf (x)[a-p + f3m + U(r)]ip(x). (11) 



Here we use a potential, which we call the core-Hartree (CH) potential, defined by 

r°° dx 

10 r m{%) 



U C n(r) = a / -^-p(x), (12) 



where r m (x) =max(r, x) and 

p(r) = J2(9l(r) + fa(r)) (13) 

a 

is the sum of the charge densities of all electrons in the filled core. This potential gives 
results close to the HF potential. To illustrate the accuracy of the potential we note that 
for cesium it predicts the removal energy of the valence 6s electron to be 3.267 eV, to be 
compared with 3.466 eV for the HF potential result. (The experimental value is 3.894 eV.) 
Because we have incorporated this potential into H , the introduction of a new interaction, 
which acts like a counterterm, 

H C T = -J d 3 x^(xho^(x)U C H(r), (14) 

is required. At the level of perturbation theory we use in this paper it plays no role. 

By working with this modified Furry representation, it is possible to put atomic many- 
body calculations in a field-theoretic framework, with Feynman diagrams involving Coulomb 
photon exchange and Hqt having a direct correspondence with the formulas of MBPT [lj| . 
Physics associated with transverse photon exchange and negative-energy states enters in 
a well-defined manner. The issue of negative-energy states, which require care to prop- 
erly include in Hamiltonian formalisms, is of course correctly treated in a field-theoretical 
approach. This issue was emphasized in Ref. 6j. 

To treat the many electrons present (up to 87 for francium), we exploit the fact that 
all atoms considered here have one electron outside a filled core, denoted \Oc)- Standard 
methods from many-body perturbation theory can then be used to treat the paramagnetic 
ground states as 

\<P) = bl\O c ), (15) 



where v denotes a valence electron outside a closed-shell core. For example, cesium would 
have v being the 6s state, and |0c) a xenon-like core, with the 54 electrons filling the Is — 5s, 
2p—5p, and 3d— 4d shells. The completeness of the shells allows angular momentum identities 
to simplify the calculations, but summations over all the core states are still needed. 



III. QUARK EDM CONTRIBUTION TO d n 

We confine our attention here to the contribution of nonvanishing edms of the up and 
down quark to the neutron edm. A thorough discussion of other contributions to the overall 
edm of the neutron is given in Chapter 13 of Ref. jl|. We will present results for the lowest 
order contribution, dn \ and for the correction due to internal electric fields, dffl ■ 

A general spin-1/2 particle with edm d is described by the interaction Hamiltonian 

#cdm = i- J d 3 ^(z)a^ 75 ^(x)F^(x), (16) 

where F^ u (x) at this point includes both classical and quantized fields. We begin by taking 
the electromagnetic field to be an external constant electric field pointing in the positive z 
direction. A neutron with edm d n in such a field has a linear Stark effect, with the energy 
shift in the rest frame of the neutron when its spin is up being 

AE = -d n E ext . (17) 

To calculate the effect the edms of up and down quarks bound in a neutron have on d n , 
we need to evaluate the energy shift in the same external electric field using the static well 
model described above, and we will then identify the coefficient of —E eXt as the neutron 
edm. Because the neutron has both up and down quarks, we generalize Eq. [16] to 



•di{x)cTQ^ip di (x) (18) 

where we have replaced F^ v with its external-field value. It is then straightforward to apply 
the Gell-Mann-Low formalism to calculate the energy shift and identify the lowest-order 
contribution as 

<#> = 0.8265 [~d d -~d u \. (19) 

The factor 0.8265 is the integral 

1 = jyr[ g l{r)+ l -f g {r)} (20) 
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for the spherical well model with massless quarks. The same kind of integral, which apart 
from the factor of 1/3 is the normalization integral, is present also in the atomic edm problem 
for valence s\/2 states, but it is extremely close to unity in that case, even for heavy atoms. 
Here, however, because the quarks are totally relativistic, the integral changes by more than 
17 percent from unity. We note that in nonrelativistic quark models the factor is taken to be 
exactly one. This 17 percent change from the nonrelativistic limit will turn out to dominate 
by far the effect of the internal electric fields in the neutron. 

We now turn to the evaluation of these effects, which are associated with one photon 
exchange. The motivation, as mentioned above, for evaluating the effect of such terms 
comes from the case of atomic physics. As will be shown in the next section, in heavy 
atoms these terms produce large enhancement factors. Because this latter enhancement is 
a relativistic effect, it is difficult to form an intuitive understanding of it. We note that 



while several explanations have been given, a recent paper [12] raises questions about their 
validity, and gives a quite different derivation involving Lorentz-Fitzgerald contraction. It 
is our view that in such situations there is no substitute for an explicit calculation of the 
effect. 

The electromagnetic field is now split into two parts, the classical field used in the lowest- 
order calculation described above, and a quantized field. This leads to three additional 
Hamiltonians, 

Hn = Qu^2 J d3x ^i{x) r r fl ^ui{x)A t '(x) +Qd^2 J d 3 x ip di (x)'j ll ip di (x)A> 1 (x) 

His = ~ ^ j ^ x [ d n^ui{x)(T a ^^ ui (x) + d d ^ di (x)a a pj 5 ijj d i(x)]F al3 (x). (21) 

In the above A 1 and are now understood to be quantized fields. The evaluation of the 
diagrams we consider, shown in Fig. [TJ requires a treatment of the quark propagator, which 
when spherical symmetry is present can be written as a partial-wave expansion in terms of 
k. This can be done with simp le modifications of a numerical approach developed for atomic 



physics described in Ref. 13] . In that paper a finite basis set using piecewise polynomials 
was applied to solving the Dirac equation for atoms. To discretize the spectrum the atoms 
were taken to be in a well with radius much larger than the size of the atom, with MIT bag 
boundary conditions applied at that radius. A typical application for atoms generates 50 
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positive-energy and 50 negative-energy states for each value of k, with the advantage that 
the first few positive energy-states accurately reproduce the known positive-energy bound 
states. All that is needed to adapt the code to the present problem is to change units and 
to eliminate the potential, as there is now no nucleus, and the screening from the other 
quarks is treated perturbatively. We emphasize that in this approach the actual excited 
states of the neutron play no role: we are simply evaluating Feynman diagrams in a simple 
bound state treatment of the neutron, using an accurate numerical method to evaluate 
quark propagators. This point will be returned to below when we note that in some terms 
the k — 1 part of the quark propagator vanishes when negative- and positive-energy states 
cancel each other, an effect which would be missed if we worked with only positive-energy 
states. 



Fig. [2] shows a part of Fig. 1(a) that involves d u and two factors of qa, with an energy 
shift we call AE a i. We discuss its evaluation in detail to illustrate how both the neutron 
and the atomic calculations are carried out. The photon and quark propagators are treated 
as follows. The photon part of the diagrams arising from a contraction of Hn and if/5 is 

(0\T(A^x)F^(y))\0), (22) 

which can be obtained from the usual Feynman gauge propagator 

/d A k e~ ik '( x ~ y ^ 
(27r)4 k2 + ig C 23 ) 




(a) (b) 

FIG. 1. The two classes of one-photon-exchange diagrams contributing to the edm of the neutron. 
The triangle vertex represents the external electric field, while the crossed vertex represents the 
edm interaction. 
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by differentiation. Note that the time component of the derivative operator acts on the 
theta functions in the time-ordering operator, but because the equal-time commutator of 
photon fields vanishes one has 

J^Xi^K - 9^ f k2+ . § . (24) 

The quark propagator in any time-independent external field, here the static well, can be 
expressed as 

(0\Tty(x)$(y))\0) = iS F (x,y) 



i 



, . dE c -iE( X0 -y ) V- ^ m {x)^m{y) , s 

2tt ^ E - e m ± iS K } 

with a factor iS for positive-energy states and — i5 for negative-energy states. The infinite 
sum over both positive- and negative-energy states is approximated as a finite sum using 
the finite basis set techniques described above. 

Using these forms for the propagators Fig. 2 corresponds to the energy shift 

/f (ftfc l e ik-(x-z) 
d 3 xd 3 yd 3 z / ? -— V k a 

(n\ : ^(f)7o7/3V'm(x)V'm(y)l/3^(f)^l(?)7oO ra/3 75^(^) : \n). (26) 

At this point we make the approximation of setting = 0, which corresponds to accounting 
for internal Coulomb fields only and ignoring internal magnetic fields. We note that in the 
atomic case, Ref. [6J did account for magnetic terms and showed they were extremely small. 




FIG. 2. Diagram QJa) for a particular assignment of quark flavors. 
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For the neutron we will find very small results for the Coulomb term, and while in this more 
relativistic case the magnetic terms may not be as suppressed as they are in a heavy atom, 
they should also be very small. The integral over d 3 k then can be written as a gradient of 
the standard integral that gives l/\x — z\. This was dealt with in two ways. In the first, a 
partial integration was carried out, which gives 

AE = 2qjd u E cxt r d 3 xd 3 yd 3 z 1 
" ' Air I \x — z\ ^-^ e r . 



(n| : ^(x^ m (x)^ m (y)y3My)^z ■ (^ u (z) l0 t^ u (z)) ■ \n). (27) 
A useful identify for evaluating the divergence term is 

V • i) m (x)E^ n (x) = ^2 IjM(-KnTn n , n m m m )Y JM (x)S mn (x), (28) 

X JM 

with 

Smn(x) = -{g m {x)g n (x))' + {fm{x)f n (x))' - — + Kn (g m (x) g n (x) + f m (x)f n (x)) (29) 

X 

and 

IjM{Kiiii,Kjfj,j) = J dVtYj M {^)x ] {^i)x{^j^j)- (30) 

The second evaluation did not use the partial integration, but applied the gradient operator 
directly to the l/\x — z\ term. These two approaches gave the same numerical final result, 
but it should be noted that in general care is needed with partial integration in the static 
well model. 

In evaluating the diagram, one encounters a sum over all possible magnetic quantum 
numbers of quarks in the initial and final neutron states, along with those of the intermediate 
quark propagator (which in fact can only have k = 1, —2 from selection rules). In the atomic 
case, these magnetic-substate sums are always complete since we sum over the filled shells 
of the core, and standard identities of Racah algebra involving products of 3j symbols can 
be used to reduce the angular portion of the calculation. Because three quarks in a neutron 
do not form filled shells, the sums in this case are not complete and such identities cannot 
be immediately applied. We therefore adopted two approaches. The first simply evaluates 
and sums the products of the three 3j symbols that enter for the various combinations 
of magnetic quantum numbers — a 'brute-force' approach. In the second, which will be 
described in detail elsewhere, the neutron wave function ffTUl) was rewritten in terms of 
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Clebsch-Gordan coefficients and a fractional parentage decomposition for the spin couplings 
of the quarks, after which the standard identities of Racah algebra could be applied after 
all. Once again, we obtained complete agreement between the two different approaches. 
Now defining 

Rij(x) = gi(x)gj(x) + fi(x)f)(x) } (31) 

and 

f R 

r ij = / dx xRij(x), (32) 
Jo 

the diagram of Fig. 2 gives 

2 K 22^ r f R f R r 

AE al = ^ad u E ext \J — r ^ E — I dx I dz^-R gm (x)S gg (z) 

e g — e m J J Q r > 

+ —ad u E ext y_. — j dx j dz^-R gm (x)S gg (z), (33) 



which leads to 

d^\al) = 1.50 x 10~ A d u . (34) 

An interesting feature of the k — 1 part of AE± a is that it vanishes, so the result above 
comes entirely from the k = —2 term. If only positive energy states are considered, a result 
on the order of the k = —2 term is found, but when negative energy states are included 
an exact cancellation takes place. The vanishing of the k — 1 term can be demonstrated 
analytically [14]. Were we to instead treat the calculation by saturating the sum with actual 
excited neutron p-states we would again obtain a nonzero result. The negative energy states 
doing the cancellation of course involve the quark sea, here treated perturbatively. 

Similar manipulations apply to the other combinations of quark lines in Fig. [IJ and the 
final formulas are, keeping only those k — 1 terms that are nonzero, 

8 Km==^-2 ^ f R f R r 

AE a = -— ot (Wd d - d u ) E ext V — ^9—1 dx I dz^-R gm (x)S gg (z) (35) 

' „ e g ~ e m Jo JO r > 

<>>m 

and 

8 K ^z^ r f R f R l 

AE b = -—a(d u - d d ) E ext } — ^— / dx dz — R gg (x)S mg (z), (36) 

which when evaluated numerically gives our result for the effect of the internal electric fields 
in the neutron, 

4 2) = 2.51 x I0~ 4 d u - 1.60 x l(T 3 d d . (37) 
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While it is of note that the contribution from the down quarks dominates that from the 
up quarks by an order of magnitude, the overall contribution from internal-field effects is 
significantly suppressed relative to the first-order neutron edm So 1 given in Eq. f ll9p . No 
effect corresponding to the dramatic violation of Schiff's theorem in heavy atoms is found, 
a result which we discuss further in the conclusions. 



IV. ATOMIC CALCULATIONS 

Turning now to the atomic case, we first note that manipulations with second-order 
perturbation theory are generally made that involve replacing the electric field felt by an 
electron in an atom, which is comprised of a one-body term coming from the nucleus and 
a sum of two-body terms coming from the other electrons, with a commutator. After 
this is done, a cancellation with first-order perturbation theory takes place, and one is left 
with an effective CP- violating one-body operator, the 'Sandars operator' [3j, given by the 
commutator 

[70S • p, d • p\ = 27075P 2 • (38) 

This commutator would vanish if the factor 70 in the first term were not present. This has 
the advantage both of showing that Schiff's theorem holds in nonrelativistic systems, when 
7o is effectively close to 1, and also of allowing one to work with only a one-body operator. 
We note that interesting discussions of alternative ways of deriving a CP-violating one- 
body operator have been given in Ref. pj, where calculations for paramagnetic atoms are 



15] , where the theory of the edm of diamagnetic atoms 



presented, and more recently in Ref. 
is treated in detail. 

If the original Sandars form is used in a first order MBPT calculation, a useful feature of 
the CH potential is that the result turns out to be exactly equal to the sum of the dominant 
three terms of the field theory calculation to be presented below. This is a special property of 
the CH potential, and use of other local potentials, while still dominated by the same three 
terms, gives results that differ from first order MBPT. Two of the three terms have already 
been encountered in the neutron calculation, but for the atomic problem an electric field 
arising from the nucleus is now present, so that we need to include an additional interaction 
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FIG. 3. A diagram representing the interaction of the valence electron's edm with the Coulomb 
field of the nucleus, represented by an X. 



Hamiltonian, 

d 3 xij(x)a a 5 ij(x)^^ — Z'{r) . (39) 

The interaction Hamiltonians used in the neutron calculation all have atomic analogs, which 
can be taken from the previous section by eliminating the sums over color, dropping terms 
involving the down quark, and replacing q u and d u with q e and d e , respectively. We retain 
the convention of keeping the valence electron spin up (the core of course has zero spin), 
but instead of presenting the coefficient of —E cxt , we give results in terms of R as defined 
in Eq. (CQ); that is, we also pull out a factor d e . We illustrate the calculation with the case 
of cesium with the CH potential, but give results for thallium and the full set of alkalis 
in Table HI The MBPT result for cesium with the standard form of the Sandars operator 
given above is an enhancement factor of 154.657. We note that this is only first-order 
perturbation theory, and the difference with more complete MBPT calculations mentioned 
in the introduction is to be expected. 

The same Feynman diagrams are present as for the neutron, but in the atomic case there 
is a sum over core states as the valence electron interacts with the electrons in the core. 
Diagrams analogous to Fig. [1] in which the photon interacts between two core states do not 
contribute to the atomic edm, because the core has total spin zero. In general, one external 
line in any Feynman diagram must correspond to a valence electron for angular symmetry 
reasons. 

The presence of the additional Hamiltonian associated with the nuclear Coulomb field 
(Eq. [39]) requires evaluation of the diagram of Fig. [3j The diagrams of Fig. [1] are now to 
be thought of as having the upper line be the valence electron and the middle line one of 
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the core states; a sum over all possible core states is implied for that line. There is also a 
second set of diagrams in which the top line is a core line and the middle line is the valence 
state. For both sets, the roles of external core and valence states on the left-hand side of 
the diagram only can also be interchanged (yielding a diagram in an 'exchange' rather than 
a 'direct' configuration). The lower "spectator quark" line is eliminated. A complication 
is the fact that there can be terms with nonvanishing energy differences flowing through 
the photon propagator (the neutron does involve photon exchange between quarks, but all 
energies are identical in that case). 

The calculations are quite similar to MBPT, with three exceptions. First, as mentioned 
above, most atomic calculations carry out a rearrangement of the Hamiltonian that builds 
in Schiff 's theorem as much as possible, leaving only a one-body operator (I38p that vanishes 
in the nonrelativistic limit, and this CP- violating operator is then treated in MBPT. Here, 
while the interaction with the nucleus is also a one-body operator, many of our diagrams are 
of two-body type, so the enhancement factor arises in a different way. Second, the electron 
propagators contain both positive- and negative-energy states. Finally, while the sums over 
positive-energy states in MBPT are generally arranged in such a way as not to include 
both occupied core states and excited states in a single sum, in our field-theoretic approach 
the electron propagator includes core states along with excited states (and negative-energy 
states). However, as long as a complete set of diagrams is included, cancellations will occur 
between different diagrams that effectively enforce this restriction where appropriate. This 
property was checked as part of the verification of our calculation. 

Because we are not building in Schiff's theorem, an important contribution arises from 
lowest order, which parallels the lowest-order neutron result. Evaluating the effect of Hj 3 
for an atom with a valence S1/2 state gives the enhancement factor 



which, as mentioned before, is very close to 1. To exhibit Schiff's theorem for a light alkali 
atom then clearly requires finding terms from the nuclear and one-photon-exchange diagrams 
that sum to a result close to —1. 

An interesting feature of the field-theoretic approach is that one automatically generates 
diagrams involving radiative corrections on either the edm or external field vertex, with a 
representative diagram shown in Fig. HI The presence of such terms was noted in Ref. jsj. 




(40) 
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Techniques to evaluate QED radiative corrections in the presence of perturbing potentials 



have been developed in Ref. 



16] , and these corrections were found to have at least a factor 



of a suppression. For that reason, we do not treat them here. 
Figure E] leads to 



R„ 



with 



dr [gi{r)gj{r) - fi(r)fj(r)] 



j 

Z(r) 



Z'(r) 



(41) 



(42) 



Since no account of the electric field coming from other electrons has yet been made, a very 
different result from the Sandars form of 154.657 is found, 



R„ 



-80.934. 



(43) 



The bulk of this comes from the 6pi/2 state, —86.354, with the next two most important 
contributions being from the 5pi/ 2 state, 7.128, and the 7pi/ 2 state, —1.182. 



Next we turn to the photon-exchange diagrams. The first diagram arises from Fig. 1(b) 
when the upper external lines are valence states and the middle external lines are core 
states. In this case, the external field interacts in the usual way with the valence electron, 
which subsequently has its edm interact with a core electron through photon exchange. It 
is this term that dominates and acts to cancel the nuclear term, as will be shown below. 



Other, smaller terms arise from Fig. 1(b) when the external core and valence lines are 



interchanged in various ways, and from Fig. 1(a) All these cases intrinsically involve two- 



body interactions, which are not particularly complicated to deal with. As with the neutron, 




(a) 

FIG. 4. Self-energy type correction to edm 
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we make the approximation of keeping only the \i = component of the photon propagator. 



The full expression for Fig. 1(b), including all combinations of external core and valence 
lines, is 



It,, = -2a J^T^ / $^tiWrA.T)t-.-U\A^r r (z)) 



ma 



+ 2a ]T^^- f (^e^-^ s -^l{x)U^ , ■ {M^M) 

J % 



ma 



+ 2a ^ -^2- / ^^e J (---)l-^l(^)^(f)V 2 • (^ m (z)£^(z)) (44) 



ma 



where v is the valence state, the states a are core states, and the sum over m is over a 
complete set of states (positive- and negative-energy) arising from the intermediate electron 
propagator. To reach this form a partial integration has been carried out, but, as with the 
neutron, we also did the calculation without doing this and obtained the same result to high 
numerical accuracy. We have dropped two terms that do not depend on v because they 
vanish by angular symmetry, as mentioned earlier, which is also the case for the second term 
in the above expression. The dominant term is the first, which can be identified with the 
charge distribution of the core interacting with the edm of the valence electron; the other 
two terms play a much smaller numerical role. 

After angular reduction, the dominant first term of R a has the form 

2<y Km= ~ Kv r r°° r°° i 

Rai = — V ^—{2j a + l) dx dz—R aa (x)S mv (z). (45) 

6 [m][a] e v~ e ™ J JO r> 

For cesium this has the value 

R al = 234.590, (46) 

again dominated by the 6px/2 intermediate state, which contributes 222.327. The accumu- 
lated sum so far is now 154.657, exactly equal to the Sandars form result. We note that 
this pattern is quite different from Ref. jo], where the nuclear term dominates. However, the 
meaning of the nuclear term is different in that work, involving a factor 70 — 1 not present 
in our approach, so the disagreement is to be expected. The other terms in Eq. (1441) . as well 



as those arising from Fig. 1(a) , are all much smaller (or zero), contributing to R at the level 



of 0.2 or less for cesium, and we will not present their values. 
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TABLE I. Atomic edm enhancement factors given by the present formalism using a CH potential 
(see text). Notation: Rq is the lowest-order term, Eq. (|4"0]) for a valence state, _R n uc is the nuclear 
term, Eq. (|4ip . and R a is the dominant photon-exchange term arising from the first term in Eq. 



Atom 
Li (2s 1/2 ) 
Na (3s 1/2 ) 
K (4s 1/2 ) 
Rb (5a 1/2 ) 
Cs (6s 1/2 ) 
Fr (7s 1/2 ) 
Tl (6pi /2 ) 



Rn 



Sum 



1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
-0.333 



-8.266 
-31.875 
-66.220 
-110.428 
-80.934 

701.730 
-573.199 



7.270 
31.314 
68.807 
143.160 
234.590 
364.161 
-219.132 



0.004 
0.439 
3.588 
33.732 
154.657 
1066.891 
-792.665 



While of course not of interest for a practical determination of d e , the case of lithium 
shows how this approach implements Schiff's theorem. The cesium results of (1, -80.934, 
234.590) go to (1, -8.266, 7.270), summing to 0.004. Three features of the lithium calculation 
are of note. The first is that the small terms we have ignored for heavy atoms, while still 
summing to a small value, can be as large as 0.7 in magnitude individually. The second is 
that we found that the approximation of dropping the exponential terms, the no-retardation 
approximation, was in some cases a poor one for a particular diagram, though this had only 
a small effect on the total answer. The third is that using potentials other than CH fails to 
exhibit the almost exact cancellation shown above for lithium. In such a case presumably 
higher-order terms would improve the cancellation, but of course atoms that do not have a 
large value of R are not of interest for putting limits on d e . In Table 1 we present the three 
dominant contributions to all alkali atoms and thallium along with their sums. We again 
emphasize that these R values correspond to first-order MBPT, and are included here to 
show the Schiff suppression for light atoms and the Sandars enhancement for heavy atoms, 
and not to represent high accuracy results. 
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V. CONCLUSIONS 



The field-theory approach given here for calculating the edm of the neutron and para- 
magnetic atoms, arising from their constituents having edms, has shown firstly that for 
the neutron, outside of a less than 20 percent normalization shift from the nonrelativistic 
result, that electromagnetic corrections are small. We attribute this smallness to the fact 
that unlike atoms, where the energy shift associated with internal electric fields is of the 
order of the energy scale of atomic transitions, in the neutron the electromagnetic energy 
shifts are of order 1 MeV, while transition energies in the static well model are hundreds 
of MeV. For the neutron, the diagrams in Fig. [I] are suppressed by the ratio of the internal 
electromagnetic interaction energy to these excitation energies [see Eqs. (156]) - (157|) ]. giving a 
result three orders of magnitude smaller than the lowest-order result (1T9"]) . Thus, although 
the quarks are highly relativistic within the model, the effect of the internal electromagnetic 
interactions on the edm is quite different from that found in heavy atoms, where the same 
class of diagrams yields large enhancement factors. 

Secondly, the field-theoretic approach has provided a different way of calculating the 
enhancement factor R, not involving the Sandars [3f operator (1551) . although both one and 
two-body effects have to be evaluated. 

We now discuss possible extensions of this approach for the case of the neutron. We 
note that most treatments of the neutron edm involve naive dimensional analysis, sum- 

n 

rule techniques or chiral perturbation theory, as described in Ref. The approach used 
here, while based on a very simple static-well model, could in principle be used to provide 
an alternative method of calculation. A particularly interesting application could be an 



(a) (b) 

FIG. 5. Representative Weinberg 3-gluon diagrams contributing to d r , 
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evaluation of the effect of Weinberg's CP-violating 3-gluon vertex 17 1 on d n . This has 



been treated by Bigi and Uraltsev 



18j and Demir, Pospelov and Ritz 19j. In the present 



approach one would evaluate the two diagrams shown in Fig. \5[ with the understanding that 
the 3-gluon vertex is the source of CP-violation and that an interaction with the external 
electric field is to be added in all lines. 

For atoms, we have demonstrated a different way of showing Schiff's theorem for nonrel- 
ativistic atoms and calculating the enhancement factor for relativistic atoms. However, the 
field theory here corresponds to only first-order MBPT, and it would be difficult to extend 
it even to second order, as ladder and crossed-ladder diagrams would have to be evaluated. 
The leading effect of these diagrams would be expected to reproduce second order MBPT, 
and arguments were given in Ref. js| that subleading terms involving negative energy states 
are numerically unimportant. Probably the most interesting extension of this work in the 
atomic case would be to radiative corrections, but because of their smallness there seems no 
need at present for the considerable effort such calculations entail. 
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